Quantum discord and its asymptotic behaviors in a time-dependent XY spin chain 
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We study the dynamics and the asymptotic behaviors of quantum discord in a one-dimensional 
XY model coupled through time-dependent nearest-neighbor interactions and in the presence of a 
time-dependent magnetic field. We find that the time evolution of the nearest-neighbor quantum 
discord in this system shows non-ergodic behaviors but is asymptotic to its steady value at the 
long-time limit. The zero-temperature asymptotic behaviors of quantum discord is only determined 
by the ratio between the coupling parameter and magnetic field, whereas the finite-temperature 
asymptotic behaviors determined by both of them. These asymptotic behaviors are sensitive not 
only to the initial values of the coupling parameter and magnetic field, but also to the final values. 
It is interesting to note that quantum discords are more robust than entanglement against the effect 
of temperature. We also find particular parameter regimes, where the nearest-neighbor quantum 
discord is enhanced significantly. 

PACS numbers: 03.65.Ta,75.10.Pq 



I. INTRODUCTION 

Entanglement [l| is a quantum correlation without 
classical counterpart and has been widely believed as the 
main reason for the computational advantage of quan- 
tum over classical algorithms. However, it has been 
found recently that there exist strong indications [2| that 
an more general quantum correlation, quantum discord 
(QD) is the resource responsible for the speed up in 

the deterministic quantum computation with one quan- 
tum bit (DQCl) 0. Moreover, quantum discord is more 
robust than the entanglement against decoherence such 
that quantum algorithms based only on quantum dis- 
cord might be more robust than those on entanglement. 
Quantum discord [1, as a up-and-coming quantum 
correlation, often arises as a consequence of coherence 
between different partitions in a quantum system, be- 
ing present even in separable states. This quantum cor- 
relation has received much attention in quantum com- 
putation 0, Q , quantum communication [3-|9| , dynam- 
ics of quantum discord [lol - [T^ . quantum phase transi- 
tions(QPTs) [Tsl - fTTt . witnessing [ll] etc. 

On the other hand, the spin chain model have been 
proposed as a very reliable model for the future quantum- 
computing technology in different solid-state systems and 
a rich model for studying the novel physics of localized 
spin systems. This spin chain can be experimentally re- 
alized as trapped ions [l^ [l^l, superconducting junc- 
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tions, coupled quantum dots [2l|, ultracold quan- 
tum gases fgll et al. The experimental progresses have 
triggered intensive theoretical research on entanglement 
and quantum discord in the one-dimensional spin chains 
[241 - 128} . Specifically, entanglement and quantum discord 
in the time-dependent spin chains is investigated [2^ 
[s^ l . The dynamics of entanglement in an XY and Ising 
spin chains has been studied by considering a constant 
nearest-neighbor coupling and in presence of a time vary- 
ing magnetic field (32i] . The entanglement dynamics in a 
time-dependent anisotropic XY model with a small num- 
ber of spins is studied numerically at zero temperature 
[33| . The time-dependent spin-spin coupling was repre- 
sented by a dc part and a sinusoidal ac part. Recently, 
Ref. [S^l carefully analyzes the time evolution of entan- 
glement in a one-dimensional spin chain in presence of a 
time dependent magnetic field h{t) and by considering a 
time dependent coupling parameter J{t). Both h{t) and 
J{t) are step functions of time [s^l- The entanglement 
undergoes a nonergodic behavior. The zero-temperature 
asymptotic behaviors of entanglement depend only on 
the ratio between the coupling parameter and magnetic 
field, whereas the finite-temperature asymptotic behav- 
iors depend on both the coupling parameter and mag- 
netic field. The aim of this work is to discuss what the 
advantages of quantum discord is in the time-dependent 
spin chain, comparing with entanglement and to explore 
the relation between the asymptotic behaviors of quan- 
tum discord and the parameter setting. 

This paper attempts to study the dynamics of quantum 
discord in the time-dependent XY model, and to explore 
the interesting asymptotic behaviors of quantum discord. 
In Sec. II, we introduce the time-dependent XY spin 



chain and describe the general solution for this model. The magnetization in the z-direction is defined as 
Sec. Ill analyzes the discord dynamics. The asymptotic 
behaviors of quantum discord are analyzed carefully in 
Sec. III. We conclude our work in Sec. IV. 

1 ^ 

II. QUANTUM DISCORD IN THE TIME j=i 
DEPENDENT XY MODEL Its expectation value is 



In this section, we briefly review the exact solution for 
the spin XY model of a one-dimensional lattices with 
N sites coupled through time-dependent couplings J{t) 
and subject to an external time-dependent magnetic field 
h{t), and introduce the geometric measure of quantum 
discord. 

The Hamiltonian for such a system is given by (33 | 
(setting h= 1) 
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N 



2 

i=l 



1=1 



(1) 

where af,a = {x,y,z} are the Pauli matrices and 7 is 
the anisotropy parameter. 

The coupling and magnetic field are represented re- 
spectively by 



J{t) - Jo + ( Ji - Jo)0{t) 
h{t)^ho + {hi^ho)e{t), 

where 9{t) is the mathematical step function 



9{t) 
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We assume that the system is initially in thermal equi- 
librium. The reduced two-spin density matrix p^^ (t) for 
this system therefore is 
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where matrix elements can be written in terms of one- 
and two-point correlation functions js^ 
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Tr[M^p{t)] 
Tr[p{t)] ' 



specifically [34 



N/2 



tanh[/3r(/io, Jo)] 



-y 



2Ji(Jo/ii - Ji/io)5pSin2[2tr(/ii, Ji)] 



+4r2(/ii, ji)(jo COS (j)p -\- ho) 



where cjjp = , 5p = 2j sin (j)p and /3 = 1/kT. k 
Boltzmann constant and T is the temperature. 

Using Wick Theorem [35], the nearest- neighbor s] 
correlation functions can be obtained as 



{S[St+i) — X Fi+i,i-\-i 
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where [34 
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Ql.m = ]y X! P COs[(to - l)(t>p] 



i{Jiho — Johi)5p sin[(m — l)<j)p] sin[4ir(/ii, Ji)] tanh[/3r(/io, Jo)] 
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Ji[Jo/ii - Ji/io]'5p sin2[2ir(/ii, Ji)] + 2r2(/ii, Ji)(Jo cos,^p + /io)| 



+ Spsm[{m ~ l)(l)p]< JoT^ihi, Ji) + 2{Jiho - Jo/ii)(Ji cos(/)p + hi) sm^[2tT{hi, Ji)] 



(12) 



Here 



III. DYNAMICS OF QUANTUM DISCORD 



T[h{t),J{t)] = {[J (t) cos ,Pp + h{t)]^ +j'^j'^{t) siii^ (t>p}^ . 

(13) 

Various measures of quantum discord and their exten- 
sions to multipartite systems have been proposed. Here 
we use the geometric measure of quantum discord [l^ , 



D{p) = mm^en„\\p ~ X\f 



(14) 



where Hq denotes the set of zero-discord states and 
|X— F| p = Tt{X—Y)'^ is the square norm in the Hilbert- 
Schmidt space. In order to calculate this quantity for an 
arbitrary two-qubit state, we write p in bloch represen- 
tation: 



I ^ ^ 3 



i=l 



(15) 

where Xi = Tr{p{ai 1)}, yi = Tr{p(U (g) Ci)} are com- 
ponents of the local Bloch vectors, Ty = Ti{p{(Ti O o-j)} 
are components of the correlation tensor, and Ui (i £ 
{1,2,3}) are the three Pauli matrices. Each state p can 
be expressed by the parameter set {x, y, T}, the geomet- 
ric measure of quantum discord is therefore given explic- 
itly M 

D{p)^]{\\x\\' + \\T\f~k,^,^), (16) 



where fc^ax is the largest eigenvalue of matrix K = xjp- 



We now come to study the dynamics of quantum dis- 
cord in the time-dependent XY spin chain. Throughout 
the paper we use N = 1000 to numerically demonstrate 
our general results, which are almost size-independent. 
We also employ three dimensionless parameters A = J//i, 
Ai = Ji/ft-i and Aq = Jo /ho for convenience. 

The XY model undergoes a second-order QPT (Ising 
transition) at the critical point (CP) Ac = 1, which sep- 
arates a ferromagnetic ordered phase from a quantum 
paramagnetic phase. When A > 1 it is claimed [13] 
there is another second-order QPT at 7c — 0, which is 
termed as the anisotropy transition. Differently from the 
Ising transition due to the external field, this transition 
is driven by the anisotropy parameter 7 and separates a 
ferromagnet ordered along the x direction from a ferro- 
magnet ordered along the y direction. 

Fig. [T] plots the dynamics of the nearest-neighbor dis- 
cord D{i,i + 1) at zero temperature for the transverse 
Ising model with 7 = 1. When the coupling parame- 
ter (and the magnetic field) is a step function, the dis- 
cord reaches a value that is neither its value at J = Jq 
(h = ho) nor at Ji (h = hi). In other words, D{i, i -|- 1) 
shows a non-ergodic behavior, which is similar to that 
of entanglement (s^ j. At higher temperature quantum 
discord remains finite a,s t — >■ 00, though the magni- 
tude of the asymptotic discord decreases. Similarly, 
D{i,i + 1) also shows the non-ergodic behavior in the 
partially anisotropic XY model with 7 = 0.5, though the 
equilibrium time is longer than that in the Ising model. 
On the contrary, D{i,i + 1) in the isotropic XY model 
(7 = 0) becomes a constant determined by the values of 
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FIG. 1. Dynamics of the nearest-neighbor discord D{i,i + 1) 
with 7 = 1, kT — and (a) ho — hi = 1 for various values of 
Jo and Ji ; (b) Jo = Ji = 1 for various values of ho and hi . 



and T = 0. We first consider the case when the system 
is initially prepared in a state with Aq = Ac, i.e., Jp = 1, 
in Fig. 2(a) When Ai — 0, the discord oscillates in time. 
The magnitude of the discord increases with Ai until it 
reaches its maximum at Ac. Interestingly, when Ai ex- 
ceeds Ac, quantum discord behaves differently from that 
of entanglement shown in [33] . The discord decreases and 
becomes steady when J dominates over h, while entan- 
glement does not. On the other hand, when the system 
initially is in the parameter region Aq > Ac, the max- 
imum of quantum discord is much smaller as shown in 
Fig. 2(b) where Jq = 5. We see that there are two peaks, 
and the second is higher and decreases with increase in 
Ao . For the partially anisotropic cases, quantum dis- 
cord behaves similar to that in the Ising case but with 
smaller magnitudes. Quantum discord of the completely 
isotropic XY system is independent of Ai. 



IV. ASYMPTOTIC BEHAVIORS OF 

Jo and Hq. The isotropy of the initial coupling param- QUANTUM DISCORD 

eters may make spins equally aligned into the x and y 
directions, apart from those in the z-direction, which re- 
sults in finite quantum discord. Increasing the coupling 
parameters strength would not change the associated dis- 



cord. 




X 




FIG. 2. (Color online)Z)(i, i -I- 1) as a function of Ai and t at 
kT = with 7 = 1, /lo = fti = 1 and (a) Jo = 1; (b) Jo = 5. 

We now focus on the dynamics of quantum discord as 
a function of Ai. In Fig. [21 we plot D{i, i -I- 1) as a func- 
tion of t and Ai for the Ising model where ho = hi = 1 



FIG. 3. D{i, i + 1) as functions of A for h — ho = hi and 
J = Jo = Ji at (a) fcr = with different J and h; (b) kT = 1 
with ho = hi ^ 0.25, 1, 4; {c)kT = 1 with Jo = Ji = 0.25, 1, 4; 
(d) kT ^3 with ho ^ hi = 0.25, 1, 4 with 7 = 1. 

We will examine the asymptotic behaviors of the quan- 
tum discord D{i, i + as a function of X a,t {t ^ oo) for 
different values of J = Jq = Ji and h = ho = hi and at 
different temperatures. First, for the Ising model (7 = 1), 
Fig. [Sja) shows that the zero-temperature behavior of 
D{i,i + 1) depends only on the ratio J/h {— A) rather 
than their individual values. D{i,i + 1) starts at zero, 
reaches a maximum and then vanishes at larger values of 
A. The maximum of asymptotic value of D{i,i + 1) de- 
creases as the temperature increases as shown in Fig.|3{b) 
and(d). On the other hand, the finite-temperature dis- 
cord D(i,i + I) is not only determined by the ratio of J 
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and h but also individual values of J and h. Fig. |3l^b) 
and Efc) indicate that an increase in h and J causes the 
maximal discord to increase when kT = 1. 
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FIG. 4. D{i, i + 1) as a function of A for h = ho — hi and 
J = Jo = Ji at (a) kT — with different combinations of J 
and h; (b) kT = 1 with ho = hi ^ 0.25, 1,4; (c) kT = 1 with 
Jo = Ji = 0.25, 1, 4; (d)fcT = 3 with /io = fei = 0.25, 1, 4 with 
7 = 0.5. 



Next, we discuss the partially anisotropic XY model 
with 7 = 0.5. Similarly to the Ising case, the zero- 
temperature discord D(i,i + 1) depends only on the ratio 
J/h as in Fig. IH^a). D{i, i + 1) starts from zero, reaches 
a maximal value and then decays to a constant value 
for larger A's, where entanglement vanishes [s^]- When 
h ^ J, the magnetic field dominates such that spins are 
aligned along the z direction and the nearest-neighbor 
discord is zero. When h <^ J, the coupling will dominate. 
For the partially anisotropic model, the strong nearest- 
neighbor couplings J make the spins aligned isotropically 
and consequently the discord maintains an equilibrium 
and finite value ( see also the results for entanglement 
[34}). The critical behavior of quantum discord around 
the critical point A = 1 changes considerably as the tem- 
perature and other parameters. The maximal discord is 
enhanced at high magnetic fields and the stronger cou- 
pling J. The temperature effect is shown in Fig. |H[d). 
The critical behavior of quantum discord disappears for 
some values of h and J. It manifests that the thermal 
excitations suppress quantum effect and even ruin the 
critical behaviors of quantum discord. 

Likewise, D{i, i + in the isotropic XY model depends 
only on the ratio J/h at T = 0, as shown in Fig. EJa). 
Differently, D(i,i -f 1) starts from zero and saturates at 
A = 8. There seems to be a tiny peak at A = 1 before 
the saturation. Interestingly, raising the temperature de- 
lays the saturations of the discord but does not reduce 
their amplitudes, as shown in Fig. [5jb). The magnetic 
fields individually affect the saturations. The smaller the 
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FIG. 5. D(i,i 4- 1) as a function of A for h = ho — hi and 
J = Jo = Ji at (a)fcr — with different combination of J 
and h; {h)kT = 1 with ho = hi ^ 0.25, l,4;(c)A:r = 1 with 
Jo = Ji = 0.25, 1, 4; (d)fcr = 3 with ho = hi = 0.25, 1, 4 with 
7 = 0. 



magnetic fields are, the later the saturations appear. In 
Fig. [5{c) , the peaks at critical point become visible when 
increasing J and noticeably the saturations oi D[i,i~\-1) 
happen with larger amplitudes. The higher tempera- 
tures delay the saturations of the quantum discord as 
in Fig.ElJd). 

To further study the roles of the magnetic field and 
coupling parameters /iq, hi, Jq and Ji, we first calculate 
the zero-temparature asymptotic behaviors of quantum 
discord at i 00 in wider parameter regimes. Fig. 12] has 
four discord contours for the Ising model at t — c». We 
plot the discord contour function of Jo and Ji in Fig.jBJa), 
where ho = hi ^ 1. The discord starts with a zero value 
for Jo = Ji = 0, reaches the maximum at Jo ~ Ji « f , 
and vanishes in the parameter regimes Jo ^ 2 and Ji > 4. 
It is interesting to note that for some values in the region 
Jo < 1 and Ji > 1, the asymptotic discord has finite 
values whereas entanglement vanishes |34l | for all values 
in this region. Fig. inijb) is the discord contour of ho and 
hi with Jo = Ji = f. The discord starts with zero at 
hi = Hq = and reaches the maximum at /lo ~ /ii ~ f . 
We then plot the contour D{i,i -I- 1) as a function of Jo 
and ho in Fig. ^c) . The discord has its maximum when 
Jo ho and disappears when Jo deviates largely from ho- 
We also show the asymptotic behavior of Z?(i, i -I- f ) as a 
function of Ji and hi as shown in Fig.[6l[d). The largest 
discord is reached at Ji = /ii. The profile of the discord 
contour versus Jo and ho is smooth, whereas the profile 
versus Ji and hi is not. The reason is that Jo and ho 
only affect the system at the initial time, whereas Ji and 
hi affect the system all the time. A little change in the 
values of Ji and hi has a great impact on the asymptotic 
behaviors of quantum discord. Fig.[7]is the same as Fig. [5] 
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FIG. 6. (Color online)The asymptotic behavior of D{i,i + 1) as a function of (a) Jo and Ji with ho — hi = l,(b)fto and hi 
with Jo = Ji — 1, (c) ho and Jo with hi — Ji = 1 and (d) hi and Ji with /lo = Jo = 1 at kT — with 7 = 1. 



except 7 — 0.5. The asymptotic behavior of discord is 
similar to that of the Ising case as shown in Fig. [71 The 
differences are that the quantum discord does not vanish 
in the range of Jq > 1 and Ji > 1, and the maxima of 
quantum discord in the partially anisotropic XY model 
are smaller than those in the Ising model. There is no 
peak in the vicinity of Ji — hi. Fig.[5]is for the isotropic 
XY model, where D{i, i + depends only on the Jq and 
ho but not on Ji and hi. 

The temperature effect on quantum discord is shown 
in Fig. [9] In Fig. IH^a), the quantum discord reaches the 
maximum at A = 1. As the temperature increases or 
A diverges from the critical value, the quantum discord 
decays more slowly than entanglement [SJ]- The asymp- 
totic behaviors of Z?(i, i + 1) as a function of Ai, Aq and 
kT, are depicted in Fig. Mjo) and IH^c). These figures 
demonstrate that the quantum discord is more robust 
than entanglement against the effect of temperature. 

Finally, we explore the asymptotic behaviors of quan- 
tum discord in the A1-7 phase space for the different 
values of Jq. In Figs. [TUlJa), [TUTb) and \W[ c) . we plot 
D{i,i 4- 1) as a contour function of Ai and 7 for given 
magnetic fields ho — hi — 1 at kT = 0. When Jo < 1, 
there is no bigger enhancement in D{i,i + 1). When Jq 
is in the vicinity of 1, there are two peaks. For the larger 
Jo, two peaks merge. When Jo > 1 and Ai < 1, D{i, 
is enhanced significantly in the 7 = region, whereas for 



Ai > 1, D{i, i + is not enhanced. 



CONCLUSIONS 



We have investigated the dynamics and the asymp- 
totic behaviors of quantum discord in one dimensional 
XY model coupled through a time-dependent nearest- 
neighbor coupling and in the presence of a time- 
dependent magnetic field. The system shows non-ergodic 
and critical behaviors. The zero-temperature asymptotic 
behaviors of the system at the long time limit depends 
only on the ratio between the coupling and the mag- 
netic field but not their individual values. On the con- 
trary, these behaviors do rely on individual values of the 
coupling and the magnetic field at the finite tempera- 
ture. The asymptotic behaviors are sensitive not only to 
the initial values of the coupling and the magnetic field, 
but also to the final values of the coupling and the mag- 
netic field. We have demonstrated that quantum discord 
is more robust than entanglement against the effect of 
temperature. This robustness will be useful in design of 
fault-tolerant quantum algorithms and in other quantum 
information processing. 
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FIG. 7. (Color onlme)The asymptotic behavior of D(i,i + 1) as a function of (a) Jo and Ji with ho — hi — 1, {h)ho and hi 
with Jo — Ji — 1, (c) ho and Jo with hi — Ji = 1 and (d) hi and Ji with ho — Jo — 1 at kT — with 7 = 0.5. 
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